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00 ■ ABSTRACT 

In Paper I, we followed the evolution of binary stars as they orbited near the supermassive black 
hole (SMBH) at the Galactic center, noting the cases in which the two stars would come close enough 
together to collide. In this paper we replace the point-mass stars by fluid realizations, and use a 
smoothed-particle hydrodynamics (SPH) code to follow the close interactions. We model the binary 
components as main-sequence stars with initial masses of 1, 3 and 6 Solar masses, and with chemical 
composition profiles taken from stellar evolution codes. Outcomes of the close interactions include 
O i mergers, collisions that leave both stars intact, and ejection of one star at high velocity accompanied 

by capture of the other star into a tight orbit around the SMBH. For the first time, we follow the 
evolution of the collision products for many (> 100) orbits around the SMBH. Stars that are initially 
> . too small to be tidally disrupted by the SMBH can be puffed up by close encounters or collisions, with 

the result that tidal stripping occurs in subsequent periapse passages. In these cases, mass loss occurs 
■ ■ ' episodically, sometimes for hundreds of orbits before the star is completely disrupted. Repeated tidal 

flares, of either increasing or decreasing intensity, are a predicted consequence. In collisions involving 
CO ■ a low-mass and a high-mass star, the merger product acquires a high core hydrogen abundance from 

£> ' the smaller star, effectively resetting the nuclear evolution "clock" to a younger age. Elements like Li, 

ON , Be and B that can exist only in the outermost envelope of a star are severely depleted due to envelope 

■ ejection during collisions and due to tidal forces from the SMBH. Tidal spin-up can occur due to 
either a collision or tidal torque by the SMBH at periapsis. However, in the absence of collisions, tidal 
spin- up of stars is only important in a narrow range of periapse distances, r t /2 < r per < r t with r t the 

OO ■ tidal disruption radius. We discuss the implications of these results for the formation of the S-stars 

■ and the hypervelocity stars. @ 

Subject headings: black hole physics-Galaxy:center-Galaxy:kinematics and dynamic 

* * 

; 1. INTRODUCTION S-stars f Perets et al. 1 120091 : iPerets fc Gualandris 1120101 ). 

X : Tidal breakup of binary stars by the supermas- Th ^ ovi ^ of the . 1 l ?™ ar 7 Progenitors of the HVSs is 

sive black hole (SMBH) at the Galactic center (GC) not clear. One possibility is that the binaries originated 

has been invoked to explain a number of otherwise at distances of a few parsecs from the GC and were 

puzzling discoveries, including the hypervelocity stars subsequently scattered inward by massive-perturbers 

HVSs) that a re observed in the halo of the Milky Way ( |Perets et al. | I2QQ2D- In this scenario, most of the bi- 

.■BrownetaUl200a [2(M 120071 iSh . and the S-stars, nanes will he on either unbound or weakly-bound or- 

apparently young, main-seque nce stars in tight eccen- blts ™ th respect to the SMBH, and they wi l encounter 

trie orbits arou nd the SMBH (lEisenahauer et al. 2005; J* only once before being scattered onto different or- 

IGillessen et, al.ll2T)0l . As first pointed out by J. Hills, bits Alternatively, the binaries may form closer to the 

close passage of a binary star near a SMBH can result in S , MB H> P erha PS ™ the young (or an older) stellar disk 

an exchange interaction, such that one component of the thaMs observgd_between - 0.04 P^d_^5_p^om 
binary is ejected with greater than escape velocity while 
the other star is scattered onto a tightly-bound orbit 
(Hills"! 119881 : lYu fe Tremaind 12003). The predictions of 



the SMBH dNavakshin fc Cuadra I 2005: P aumard et al. I 
20061: iNavakshin fc Cuadra 1 1200/1 : iLevinl 120071 )" Also, 



Perets I (|2009bl ) suggested that binaries could be left near 



this model are broadly con sistent with the observed prop- ^ e G , C trough a triple disruption by the SMBH. In 
erties of both the HVSs (iBromlev et a l. 2006) and the tti^e latter cases the binaries would be bound to the 

SMBH and would encounter it many times before being 
i — . — — . disrupted. 

I ant onmi O astro . r it . edu I Tr , T n •> r , , i ■ , , 

- Movies can be found at htt p://a S trophysic S .rit.edu/fantonini/tbbs2/| If the fimte slz es of stars are taken into account, a 
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number of outcomes are possible in addition to simple 
binary disruption. The two stars can collide, resulting in 
a merger i f the relative velocity is l ess than stellar escape 
velocities (jGinsburg fc Loebl 12007). Since the radius of 
tidal disruption of single stars by the SMBH is compa- 
rable to the binary disruption radius, stars can also be 
tidally disrupted by the SMBH, either before or after 
their close int eraction with each ot her. 

In Paper I (j Antonini et al.|[2010f ) we presented the re- 
sults of a large number of Af-body integrations of point- 
mass binary stars on eccentric orbits around the GC 
SMBH. In many cases, the trajectories of the two stars 
were found to imply a physical collision, assuming that 
the unperturbed stars had radii similar to those of nor- 
mal main-sequence stars of the same mass. The proba- 
bility of physical encounters was found to increase sig- 
nificantly if the binaries were allowed to complete many 
orbits about the SMBH. In some cases, one or both stars 
also passed close enough to the SMBH that gravitational 
tides would be expected to significantly affect their in- 
ternal structure. 

In this paper, we use smoothed-particle hydrodynam- 
ics (SPH) simulations to study the binaries from Paper 
I that approached closely enough to physically interact. 
The A'-body simulations of Paper I were first used to 
identify initial conditions that resulted in close interac- 
tions between the two stars. The point-mass stars were 
then realized as macroscopic, fluid-dynamical models and 
integrated forward in the gravitational field of the SMBH 
using an SPH algorithm. As in Paper I, we followed the 
trajectories for multiple orbits around the SMBH, allow- 
ing us, for the first time, to investigate the consequences 
of repeated tidal interactions with the SMBH. 

In §2 we briefly discuss time scales for binary disrup- 
tion at the GC. Our initial conditions and numerical 
methods are described in §3 and the results in §4. Some 
observable consequences are presented in §5. §6 sums up. 

2. THE SURVIVAL TIME OF BINARIES AT THE 
GALACTIC CENTER 

In a dense environment, binaries may evaporate due to 
dynamical interactions with field stars if 



\E\/(M h a 2 ) < 1, 



(1) 



with E the internal orbital energy of the binary, Mb 
the binary mass, and a the one-dimensional velocity dis- 
persion of the stellar background. In principle, because 
most of the binaries at the GC are expected to be "soft" , 
\E\ < Mb<r 2 , and because the binary evaporation time 
t ev is a function of the distance from the SMBH, the 
variation of t ev with galactocentric r adius can be u sed 
to constrain the origin of the HVSs (jPerets I l2QQ9ah . If 
the evaporation time at some radius is shorter than the 
lifetime of a typical main-sequence star, the stellar popu- 
lation in this region would be dominated by isolated (i.e. 
single) stars. 

Here we show that the survival time of binaries at 
galactocentric distances r < O.lpc is likely to be com- 
parable to the typical main-sequence lifetimes of most 
stars in this region. We also show that, within a radius 
of r ~ 0.3pc, t ev becomes essentially independent of ra- 
dius. 

Beyond ^ 1 pc from Sgr A*, the mass density de- 
termined from the stellar kinematics follows p ~ r' 




0.01 0.1 



Radius(pc) 

Fig. 1. — Evaporation time of binaries vs. galactocentric radius 
for different values of the binary semimajor-axis ao- Solid curves 
show the evaporation time for the density model of equation (0 
with 7 = 0.5 while the dashed curves correspond to the coreless 
model with slope 7 = 1.8. The filled grey region gives the ages of 
the S-stars (Eisenahauer et al. 2005). 



1-5 < /3 < 2 (e.g. IQh et al. I l2QQ9h , At smaller radii, 
number counts of the dominant (old) stellar population 
near the GC suggest a space density that is weakly 
rising, or falling, toward the SMBH, i nside a core of 
radius ~ 0.5pc (Ifi uch holz et al.l 120091 : IDo et al.l 120091 : 
iBartko et al.ll2010f ). Approximating the mass density as 
a broken power-law, 



p( r ) = Po ( - 



1+1^ 



(7-/3)/2 



(2) 



with ro = 0.3pc and (3 = 1.8, and setting p = 1-3 x 
lO 6 M0pc -3 g ives a good fit to the space density outside 
the core (e.g. iMerritt Il2010f ). At smaller radii, the un- 
certainties in p are represented by the poorly determined 
value of 7. 

Th e evaporation time is given by (jBinney fc Tremaind 



M 16v/7r/9a lnA' 



(3) 



-/3 



where InA is the Coulomb logarithm, M the mass of the 
field stars, ao the binary semimajor-axis, and a is calcu- 
lated from the Jeans equation, 

poo 

p(r)a(r) 2 =G dr'r'' 2 [M. + M*(< r')} p(r'), (4) 

J r 

with M*(< r) the total mass in stars within r, and M. 
the mass of the centr al black hole . Her e after, we adopt 
M. = 4 x 10 6 M (jGhez et al. I [20081 : IGillessen et "all 
2009). In Figured] we plot the evaporation time of bina- 
ries in the density model of equation fl5J) as a function of 
galactocentric radius, assuming Mb = 2M, In A = 15 and 
two different values of the internal slope: 7 = 0.5 repre- 
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Fig. 2. — Fractional chemical abundances (by mass) versus en- 
closed mass fraction m/M for our M = IMq (red curves), 3M0 
(green curves), and 6Mq (blue curves) stars, as calculated by the 
TWIN stellar evolution code. 
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Fig. 3. — Internal energy U, gravitational potential energy W, 
kinetic energy T, and total energy E versus time t for the relax- 
ation (left panels) and subsequent dynamical evolution in isolation 
(right panels) of the SPH model for a 6 Mq star. Note that the 
time t is shown on different linear scales for the relaxation and the 
dynamical evolution. Energies are in units of 10 48 erg. 



sentative of the observed distribution and 7 = 1.8 which 
correspo nds approximately to a relaxed system around a 
SMBH (iBahcall fc Wolfl[l976h . 

From the figure it is clear that the survival time of 
binaries would be greatly increased if the distribution of 




12 3 
r 



Fig. 4. — Radial profiles of the SPH model for a 6 Mq star both at 
the end of relaxation (upper panels) and after 4200 days of hydro- 
dynamical evolution (lower panels). The frames in the left column 
show profiles of pressure P, density p, temperature T (in Kelvin), 
and mean molecular weight \i in units of the proton mass ra p , with 
the dashed curve representing results the TWIN evolution code and 
dots representing particle data from our SPH model. The right col- 
umn provides additional SPH particle data: individual SPH parti- 
cle mass m-i, smoothing length hi, number of neighbors Njy, and 
radial component of the hydrodynamic acceleration ahydro (upper 
data) and gravitational acceleration g (lower data). Unless other- 
wise stated, quantities are in solar units (G = Mq = Hq = 1). 



stars at the GC has a low density core (for comparison, 
see also Figure 1 in Perets [2009]). The figure also shows 
that, at any radius, for ao < 1 au, this time is larger 
than the typical lifetime of the the S-stars. Therefore, it 
cannot be excluded that the S-stars were initially part of 
binary systems originating at galactocentric distances of 
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few tens of milliparsecs. 

We finally note that, in the context of this paper, 
it might be more appropriate to compare t ev with the 
time scale required to drive the eccentricities of the 
stars away from their initial values to produce quasi- 
radial orbits. At the GC this time can be of the or- 



der of few Myr ( Lockmann et al. 2008; Madigan et al. 



2008; Merritt et al. 2009; Perets et al. 2009; Fujii et al 



20101 ). which is much shorter than the typical evapora- 
tion time scale of binaries with ao < 0.2 au. 

3. INITIAL CONDITIONS AND NUMERICAL 
METHOD 

3.1. Orbital initial conditions 

In Paper I we used the high- accuracy numer ical inte- 
grator ARCHAIN (jMikkola fc Merritt l2QQ8L l2QQ6) to study 
the dynamics of main-sequence binary stars on highly el- 
liptical orbits whose periapsides lay close to the SMBH. 
We determined the final orbital properties of both ejected 
and bound stars. Initial conditions consisted of equal- 
mass binaries on circular relative orbits with random 
orientations, initial separations ao in the range 0.05 - 
0.2 au, and individual stellar masses M of 3 M — 6 M . 
The binaries were given a tangential initial velocity in 
the range 4 km s _1 - 85 km s _1 and initial distances of 
d = 0.01 — 0.1 pc from the SMBH. Using the mass-radius 
relation R/R Q = (M/M ) 0,75 ([Hansen et all 12004( 1. we 
could assign a physical dimension to the particles and 
investigate the probability of stellar collisions and merg- 
ers. In detail, we defined the minimum impact parameter 
for a collision as 2i?, and the two stars were assumed to 
coalesce when their relative velocity upon collision was 
lower than the escape velocity from their surface. The 
binary separations adopted in this work were close to 
the extremes of the interval within which HVSs can be 
produced: small semimajor-axes, ao < 0.05 au, result in 
contact binaries, while for ao > 0.2 au few stars would 
be ejected with velocitie s sufficient to escape the Galaxy 
(jGualandris et al.ll2005f ). 

In order to treat also the case of unequal- mass binaries, 
we extended the work of Paper I to include an additional 
set of ~ 2000 integrations of binaries with component 
masses M2 = 1 — 3M and Mi = 6M . Initial condi- 
tions and results of these new runs are summarized in 
Appendix A. In the following, for the case of unequal 
mass binaries, we distinguish between the primary and 
secondary stars' quantities using the subscripts 1 and 2 
respectively. 

The orbital initial conditions that we adopt in the SPH 
simulations correspond to binaries that enter well within 
their tidal disruption radius rbt, approximated by the 
expression (jMiller et al. 1120051 ): 



M. 

Mb 



1/3 



ao, 



(5) 



with Mb = Mi + M2 the total mass of the binary. There- 
fore, the binaries in all of our simulations are strongly 
perturbed at the first periapse passage. Moreover, we 
focus on cases where the tidal perturbations from the 
SMBH on the single stars are expected to be signifi- 
cant. This corresponds to binaries with periapsides that 



lie close to the tidal disruption radius of a single star, or 

(6) 



n 



We selected three types of initial condition that can 
be classified according to the final outcome of the SPH 
simulations. 

1. Stellar collision (without merger). In some cases, 
gravitational perturbations from the SMBH can lead to 
a physical collision between the two stars. If the relative 
velocity at collision is sufficiently high, the stars can sur- 
vive the interaction and avoid a merger. Soon after the 
collision, one star can be ejected at high velocity. 

2. Stellar merger. If the two stars collide with a relative 
velocity at impact smaller than approximately the escape 
velocity from their surfaces, coalescence occurs, resulting 
in the formation of a new, more massive star. The merger 
remnant will remain bound to the SMBH unless extreme 
mass loss occurs during coalescence. 

3. (Clean) ejection of a HVS. The tidal breakup of the 
binary by the SMBH results in the ejection of a star at 
very high velocity. The former companion of the ejected 
star loses energy in the process and is deposited onto a 
tight orbit around the SMBH. Here, "clean" means that 
the member stars of the binary do not collide with each 
other during the process of ejection. 

For the sake of simplicity, in all the SPH simulations 
we chose the initial apoapsis of the external binary orbit 
to be d = 2000 au ~ 0.01 pc (except for case H8 which 
has d = 1700 au; see Table 1). However, starting the SPH 
simulations from this distance would greatly increase the 
required computational time. The stars were therefore 
initially placed at a point of the orbit corresponding to a 
much smaller distance from the SMBH: ro = 5r*bt- This 
choice for ro was motivated by the fact that at these 
distances, the tidal forces from the SMBH are still too 
weak to significantly influence the internal structure of 
the stars. In addition, since ro is considerably larger than 
rbt, at the initial time the internal binary eccentricity is 
still close to zero (e < 10 -2 ). As in Paper I, the plane of 
the binary's orbit about the SMBH is the x — z plane. 

3.2. SPH numerical techniques and initial conditions 

SPH is a Lagrangian method in which the fluid is rep- 
resented by a finite number of fluid elements or "parti- 
cles." Associated with each particle i are, for example, 
its position r*, velocity v*, and mass ra$. Each parti- 
cle also carries a purely numerical smoothing length hi 
that determines the local spatial resolution and is used 
in the calculation of fluid properties such as ac celeration 
and d ensity. For a recent review of SPH, see Ro sswod 
(2009 ). The SPH code used in this work is pre sented 
in iGaburov. Lombardi. fc Portegies Zwartl (j2010f ) , with 
the augmentation that the analytic solution to the Ke- 
pler two-body problem can be used to advance a star 
bound to the SMBH through those portions of the or- 
bit when hydrodynamic effects are negligible (see ^4.4j) . 
The equation of state is ideal gas plus radiation pressure 
and radiative cooling and heating is neglected. To cal- 
culate the gravitational accelerations and potentials, we 
use direct summation on NVIDIA graphics cards, soft- 
ening with the usual SPH kernel as in Hernquist & Katz 
(1989). Thus, gravity is softened only in interactions be- 
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TABLE 1 

Summary of the SPH simulations. 



Run 


Mi(M ) 


M 2 (M ) 


a (au) 


r per (au) 


Ai (A 2 ) 


Ci(C 2 ) 


SPH 


TV-body 


CI 


3 


3 


0.05 


1.50 


1.36 


0.668 


Collision+HVS 


Collision+HVS 


C2 


3 


3 


0.2 


2.67 


2.43 


1.19 


Collision+HVS 


Collision+HVS 


C3 


6 


6 


0.05 


1.50 


1.07 


0.526 


Collision+HVS 


Merger 


C4 


6 


6 


0.2 


2.04 


1.46 


0.716 


Collision+HVS 


Collision+HVS 


C5 


6 


3 


0.1 


5.61 


4.01 (5.09) 


1.96 (2.49) 


Collision+HVS (primary) 


Collision+HVS (primary) 


Ml 


3 


3 


0.05 


5.05 


4.59 


2.25 


Merger 


Merger 


M2 


3 


3 


0.05 


1.50 


1.36 


0.668 


Merger 


Merger 


M3 


3 


3 


0.05 


4.18 


3.79 


1.86 


Merger 


Collision+HVS 


M4 


3 


3 


0.2 


4.18 


3.79 


1.86 


Merger 


Merger 


M5 


3 


3 


0.2 


38.2 


34.7 


17.0 


Merger 


Merger 


M6 


6 


6 


0.05 


2.67 


1.91 


0.935 


Merger 


Collision+HVS 


M7 


6 


6 


0.05 


2.67 


1.91 


0.935 


Merger 


Merger 


M8 


6 


6 


0.05 


4.18 


2.99 


1.46 


Merger 


Merger 


M9 


6 


6 


0.05 


5.05 


3.62 


1.77 


Merger 


Merger 


M10 


6 


6 


0.2 


2.04 


1.46 


0.716 


Merger 


Merger 


Mil 


6 


6 


0.2 


20.4 


14.6 


7.13 


Merger 


Merger 


M12 


6 


3 


0.1 


8.08 


5.78 (7.34) 


2.83 (3.59) 


Merger 


Merger 


M13 


6 


1 


0.1 


3.48 


2.49 (5.28) 


1.22 (2.59) 


Merger 


Merger 


HI 


3 


3 


0.05 


1.50 


1.36 


0.668 


HVS 


HVS 


H2 


3 


3 


0.05 


7.07 


6.42 


3.14 


HVS 


HVS 


H3 


3 


3 


0.2 


2.04 


1.86 


0.909 


HVS 


HVS 


H4 


6 


6 


0.05 


1.50 


1.07 


0.526 


HVS 


HVS 


H5 


6 


6 


0.2 


2.04 


1.46 


0.716 


HVS 


HVS 


H6 


6 


6 


0.2 


0.60 


0.429 


0.210 


HVS 


HVS 


H7 


6 


1 


0.1 


1.44 


1.03 (2.19) 


0.503 (1.07) 


HVS (secondary) 


HVS (secondary) 


H8 


6 


1 


0.1 


0.795 


0.569 (1.21) 


0.278 (0.592) 


HVS (secondary) 


HVS (secondary) 



licity Z = 0.02. The 3M star is evolved to an age of 
50 Myr, yielding the same 2.15R© (0.01 au) radius as in 
the corresponding models of Paper I. The 1 and 6M 
stars are each evolved to 18.2 Myr, yielding stellar radii 
of 0.891 and 3.44R respectively. Figure [2] shows the re- 
sulting composition profiles for our 1, 3, and 6 M stars, 
colored red, green, and blue respectively. 

Initially, we place the SPH particles on a hexagonal 
close packed lattice, with particles extending out to a dis- 
tance only a few smoothing lengths less than the full stel- 
lar radius. After the initial particle parameters have been 
assigned according to the desired profiles from TWIN, we 
allow the SPH fluid to evolve into hydrostatic equilib- 
rium. During the relaxation calculation, the drag force 
we include is the normal artificial viscosity but in the 
acceleration equation only. 

Figure [3] shows energies versus time both during and 
after the relaxation process. From the internal energy 
U and potential energy W curves, it is apparent that 
the star oscillates on a hydrodynamical timescale, specif- 
ically with a fundamental period of about 0.08 days. 
During relaxation, these oscillations are damped by the 
drag force, which does negative work on the system and 
decreases the total energy E toward that of a mini- 
mum energy equilibrium state. At a time of 1.84 days 

(= 100G -1 / 2 Mq 1 ^ 2 Rq^ 2 ), the drag force is removed and 
the star is allowed, as a test of stability, to evolve dy- 
namically in isolation. During this dynamical evolution, 
the internal energy U and gravitational energy W each 
remain nearly constant. By t = 9 days, the kinetic en- 
ergy T has diminished to nearly 10 orders of magnitude 
less than the total energy E in magnitude, corresponding 



tween neighbors and it is softened by exactly the SPH 
kernel of each particle. As the smoothing length of a 
particle changes, so does its softening. The use of such a 
softening with finite extent (as opposed, for example, to 
Plummer softening) increases the accuracy and stability 
of SPH models, consistent with the studies of Athanas- 
soula et al. (2000) and Dehnen (2001). 

In our simulations, the SMBH is a compact object par- 
ticle that interacts gravitationally, but not hydrodynam- 
ically, with the rest of the system. The gravity of the 
SMBH is softened according to a density profile defined 
by the standard SPH cubic spline kernel with a constant 
smoothing length h m = 20R©. This approach has the ad- 
vantage that the treatment of gravity is unsoftened for 
separations r > 2h m . We note that h m is small compared 
to the periapsis separation r per in all cases, so our code is 
able to follow bound stars around the black hole for many 
orbits without introducing spurious secular effects from 
gravitational softening. The SMBH is allowed to move 
in response to gravitational pulls. However, because the 
4 x 10 6 M SMBH is much more massive than any of 
the binaries being considered, the SMBH always remains 
very near the center of mass of the system, which we take 
to be the origin. 

Before simulating the interaction of a binary with the 
SMBH, we must first prepare an SPH model for each 
binary component in isolation. To compute stellar struc- 
ture and compo sition profiles, we use the TWIN stellar 
evolution code (lEggletonl Elzl iGlebbeek fcToIsI [20081 : 
Glebbeek 20 08|) from the M USE software environment 
(|Portegies Zwart et al.ll2009[ ). We evolve main-sequence 
stars with initial helium abundance Y = 0.28 and metal- 
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to an exceedingly small amount of noise in an otherwise 
static model. The overall level of energy conservation is 
excellent: extrapolating forward the drift in total energy 
E, which is linear in time, we find it would take about 
2.4 x 10 5 days (660 years or 3 x 10 6 oscillation periods) 
of hydrodynamical evolution to reach a 1% error in total 
energy. 

Our approach allows the parent stars to be modelled 
very accurately. As an example, Figure |4] plots both de- 
sired profiles and SPH particle data for the 6M star. 
The structure and composition profiles of the SPH model 
closely follow the desired TWIN profiles. Our relaxed 
models remain static and stable when left to evolve dy- 
namically in isolation: indeed, the particle data shown 
in the lower panels of Figure |4] are nearly identical to 
those in the upper panels, demonstrating that there are 
no significant changes in the model even after more than 
4000 days (or equivalently 5 x 10 4 oscillation periods) of 
hydrodynamical evolution. 

The binaries in our simulations are then created sim- 
ply by shifting two stellar models, each taken from the 
end of a relaxation calculation, to the appropriate initial 
position and velocity provided by the N-body code. We 
begin with binary components irrotational in the inertial 
frame, which allows us to more easily study any rotation 
imparted during the subsequent interaction. 

Unless otherwise noted, our simulations employ N w 
4 x 10 4 SPH particles: such a particle number provides 
an appropriate balance between resolution and the need 
sometimes to follow the hydrodynamics for time inter- 
vals exceeding 10 5 dynamical timescales (corresponding 
to hundreds of orbits around the SMBH). 

3.3. Timescale considerations and orbital advancement 

Because of shock heating in collisions and mergers, in 
addition to tidal heating during the periapse passage, the 
bound stars are out of thermal equilibrium and larger 
than a normal main-sequence star of the same mass. 
The global thermal readjustment of the bound stars pro- 
ceeds on a thermal timescale ^thermal ~ U/L, where U 
is the total internal energy in the star and L is its lu- 
minosity. The SPH simulations confirm that the inter- 
nal energy U of the bound star after one periapsis pas- 
sage is comparable to the total internal energy of the 
star(s) from which the bound star came and typically 
U « (2 — 7) x 10 49 erg, with the larger values generally 
corresponding to more massive stars. For weak collisions 
and clean ejections of HVSs, the luminosity of the bound 
star will be comparable to the value it had in the ini- 
tial binary: the thermal timescale in such cases is then 
roughly 10 5 to 10 7 years. G uided by calculations of blue 
stragglers ([Sills et al.lll997[ ). we estimate that the lumi- 
nosity of a bound star produced by a merger or strong 
collision may be up to ~ 100 times larger than that of a 
main-sequence star of the same mass. Thus, the luminos- 
ity of our most massive merger products could be briefly 
as large as ~ 1O 6 L ^ 4x 10 39 erg s _1 , so that the global 
thermal timescale ^thermal ^ 600 years (although the lo- 
cal thermal timescale in the outer layers of the star could 
be less). We conclude that thermal adjustment over an 
orbital period is small and often completely negligible, 
and we therefore do not attempt to model the thermal 
relaxation here. 

Although the orbital period is small compared to the 



thermal timescale, it is large compared to the hydrody- 
namical timescale, which is about an hour. Following the 
full hydrodynamics of a multiple orbit encounter would 
therefore not be practical. What we do instead is wait for 
the star(s) to move sufficiently far away from the black 
hole and then advance any bound star around most of 
its Kepler two body orbit. At the same time, we remove 
from the simulation any HVS, any eject a, and any gas 
that has become bound to the SMBH. As long as the pe- 
riapse passages are treated hydrodynamically, our results 
are not sensitive to precisely which portion of the orbit is 
treated in the two body approximation. In practice, we 
wait at least 8 days after periapse, and at least 8 days 
after the merger or ionization of a binary, before mea- 
suring the orbital elements and implementing the two 
body analytic solution. The orbital advancement is per- 
formed such that the distance from the black hole to the 
bound star is unchanged but that the objects are now 
approaching one another. We preserve the orientation of 
the orbit and spin of the bound star during the orbital 
advancement. 

In order to test the reliability of the method, we run 
a portion of the first orbit for one of our simulations 
(CI in Table 1) without any orbital advancement. We 
found that, at about 20 days after we would have applied 
the advancement, their masses had decreased by an ad- 
ditional 0.002 Mq. We note that even though we retain 
too much mass, this "extra" mass is far from the stars, so 
it doesn't significantly participate in the hydrodynamics 
and it will get ejected in the next passage. We conclude 
that neglecting hydrodynamics far from periapsis is a 
reasonable approximation. 

4. RESULTS 

The first important result of the SPH simulations per- 
formed in this work is that their qualitative outcome 
agrees very well with that of the N-body integrations de- 
vised in Paper I. Among a total of 26 simulations, only in 
3 cases did the N-body approach fail to match the results 
of the hydrodynamic calculations. 

Table 1 reports the chosen initial conditions as well 
as a qualitative description of the final outcome of the 
SMBH-binary interaction in both SPH and A^-body sim- 
ulations. The strength of the SMBH-stars interaction 
is parameterized by the dimensionless quantity: A = 
^per/^t- In general, stars on orbit s meeting the con- 
dition A < 1 are tidally disrup ted (jLuminet fc Cartel"! 
1986; lEvans fc Kochanekl 119891 ). However, even when 
tidal disruption does not occur, we expect that mass will 
be stripped from the o uter regions of a ny star passing 
within its Roche limit (jPaczvnski 11197 lh . In the table, 
£ gives the Roche lobe radius (evaluated at periapsis) in 
units of the stellar radius: 



c = 



0.49<? 2/3 



O.64 2 / 3 + ln(l + ? i/3) 



R ' (7) 



with q = M/M. ([EggletoiJI 19831 ). Although this formula 
was derived under the assumption of circular orbit, it has 
been shown to work reason ably well even for e ccentric bi- 
naries, if used at periapsis (Re gos et al.ll2005l ). Note that 
because q << 1 in the present work, the approximate re- 
lation £ & 0.49 A exists between A and £. 
The results of our SPH simulations are presented in 
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what follows. We first descri be the pro duc t of o ne binary- 
SMBH interaction (Sections I4TT1 l4~2l and Q)l . and then 
we successively follow the evolution of the bound stars 
as they perform several revolutions around the SMBH 
(Section @~4]). 



4.1. Stellar Collisions 

iGinsburg fc Loebl (|2007| ) noted that the tidal breakup 
of stellar binaries interacting with the SMBH can lead, 
under some circumstances, to a physical collision be- 
tween the two member stars. In this section we inves- 
tigate the cases in which the two stars collide with a rel- 
ative impact speed large enough that they do not merge 
upon impact. Some results of the SPH simulations are 
shown in Table [2j where we also compare the asymptotic 
ejection velocity of the HVSs (i> e j), the semimajor-axis 
(a) and eccentricity (e) of the captured stars with the 
same quantities obtained in the TV-body simulations. 

The agreement between the two methods is remark- 
ably good. However, the SPH simulations systematically 
produce smaller values of v e j and larger a. This is a con- 
sequence of the efficient energy and angular momentum 
transfer occurring between the stars: the ejected star 
slows down and the captured star gains speed during the 
collision. 

From Table [2j it is clear that the mass-loss from the 
binary is smaller than a few percent of the total mass, 
and always below ~ 0.1 M . This is consistent with pre- 
vious simulations of stellar c ollisions that often found 
a small fractional m ass-loss dBenz fc Hillsl 119871 119921 : 
iFreitag fc Benzl 12005). Our calculations indicate that, 
after a collision, the mass ejected from the SMBH-stars 
system is comparable to, although smaller than, the mass 
that is ejected from the binary but remains bound to the 
SMBH; this mass-loss has a small but measurable im- 
pact on the subsequent evolution of the stars' orbits as 
demonstrated by comparing the SPH and the TV-body 
quantities in the table. Debris will eventually settle into 
a torus-like structure about the SMBH, that will sub- 
sequently evolve due to viscosity, mass inflow, radiative 
cooling, and winds. 

Spin-up is expected to be one of the main signa- 
tures of either a (off-axis) collision ([Alexander fc Kumarl 
[200lh or a tidal encounter with a massive black hole 
(Ev ans fc Kochanek|[l989l ). In our simulations the close 
stellar encounter as well as the SMBH tides at periapsis 
lead therefore to some degree of rotation in the stars with 
angular frequency: 



/G(Mi+M 2 ) GM. 



per 



(8) 



where and Q. are respectively the angular velocity 
induced by the interaction with the companion star and 
that imparted by the SMBH tides; tq is the distance of 
closets approach between the stars. The ratio in 
the cases considered here is always grater than 1 and 
varies from a maximum of ~ 20 (simulation C5) to a 
minimum of ~ 2 (simulation C4). 

Figure [5] plots the te mporal evoluti on of the dimen- 
sionless spin parameter ([Peebles Ill97l[ ): 
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Fig. 5. — Temporal evolution of the dimensionless spin parameter 
J denned in equation ([9} and the stellar masses M in solar units 
at times near the first periapsis passage in simulations with stellar 
collisions. Dashed curves correspond to the stars captured by the 
SMBH, while the dotted curves are for the ejected stars. The time 
coordinate is given in units of days and is shifted in order to have 
t = at the moment of the closest approach of the binary with 
the SMBH. The curves terminate at the time the orbit is advanced 
using the analytic two-body solution. 



where L is the spin angular momentum of the star and E 
its binding energy. In all cases, there is a sharp increase 
in the stars' spins during the periapsis passage followed 
by a second gradual decrease toward the final relaxed 
(spinning) configuration. Note that the captured stars 
have values of the final spin slightly larger than that of 
the ejected stars. This finding is consist ent with the re- 
sults of paper I (but see lSari et al.l (j201Qh as well), where 
it is shown that the captured member is always the star 
with the smallest value of the closest approach distance 
to the SMBH implying, as we expect from equation (j8]), 
a larger tidal torque at periapsis. 

Figure [5] also shows the temporal evolution of the 
stellar masses near the first periapsis passage. Each 
star typically loses about 1% of its mass, with captured 
stars (again, those that pass closer to the SMBH) losing 
slightly more mass than their ejected counterparts. In 
the cases with an equal mass binary, both stars lose a 
comparable amount of mass. The captured star in sim- 
ulation C5 actually gains mass that had been lost from 
its binary companion, as discussed in more detail below. 
In all cases, the stellar masses stabilize to an essentially 
constant value by the time the orbital advancement tech- 
nique is implemented, which is where the curves termi- 
nate. 

Figure [6] presents column density snapshots for simu- 
lation C5. The simulation models a stellar binary with 
ao = 0.1 au and with components of masses Mi = 6M 
and M2 = 3 M (see Table 1). The time indicated on the 
panels is shifted in order to have t = at the moment of 
the closest approach of the binary with the central black 
hole. Between t = and t = 0.9 days the reference frame 
is the center of mass of the binary while, in the lower- 
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TABLE 2 



Ejection velocity v e j of the HVS, orbital semimajor-axis a and eccentricity e of the captured star, and distance of closest approach between the two 
stars (ro) in the SPH (TV-body) calculations. The quantity AM^/M^ gives the fraction of mass lost from the binary, while AM 9 /M] => is the fraction 
of the mass lost from the binary that remains bound to the SMBH. All quantities are evaluated after the first periapsis passage, once the stars have 
retreated far from the SMBH. 



Run 


(km/s) 


a 

(au) 


e 


r /(R 1 + R 2 ) 


AM h /M 


b 




AM. /M b 




CI 


4608(5681) 


124(101) 


0.988(0.985) 


0.23(0.32) 


-1.33 x 


10" 


2 


7.19 x 10" 


3 


C2 


3878(3964) 


159(157) 


0.983(0.983) 


0.49(0.50) 


-2.80 x 


10" 


3 


2.06 x 10" 


3 


C3 


3335(4840) 


190(117) 


0.991(0.987) 


0.38(0.45) 


-1.24 x 


10" 


2 


6.73 x 10" 


3 


C4 


1467(1554) 


380(377) 


0.995(0.995) 


0.67(0.87) 


-1.29 x 


10" 


3 


7.08 x 10" 


4 


C5 


1764(2367) 


212(163) 


0.974(0.965) 


0.42(0.61) 


-2.16 x 


10" 


3 


1.53 x 10" 


3 




Fig. 6. — Column density plots for simulation C5 on the X — Z plane. In this case the binary has an internal semimajor-axis ao = O.lau 
and its components have masses Mi = 6Mq and M2 = 3Mq. Time t = corresponds to the periapsis passage of the binary external 
orbit. Between t = and t = 0.9 days the panels are centered on the binary center of mass. The two bottom-right panels are centered on 
the center of mass of either the captured (left) or ejected (right) star. The black hole is outside the images. At t = 0.13 days the stars 
collide. Subsequently, the 6 Mq member is ejected at hypervelocity while the secondary star remains bound to the SMBH. 



right panels, we switch to the frame in which the center 
of mass of either the captured (left) or ejected (right) 
star is at the origin. 

The first contact between the stars occurs ~ 0.1 days 
after the periapsis passage. The interaction leads to an 
episode of mass transfer between the stars (observed at 
~ 0.35 days in the figure). The smaller star gains mass 
(~ 0.02 Mq) in the collision, while the larger star loses 
~ 0.04 M . Only 6 x 10 -3 M becomes ejecta from the 
entire system (i.e., stars plus SMBH) after the first pe- 
riapsis passage, while the remainder of the gas lost from 
the binary remains bound to the black hole. The mass 
loss from the binary upon impact is therefore of order 
~ 10 -2 Mq. Note that, in this simulation, the penetra- 
tion factor A of both stars is large enough that the mass 



loss from the binary can be completely attributed to the 
stellar impact rather than to the tidal perturbations from 
the SMBH. 

The two bottom-right panels give column density plots 
of the captured (left) and ejected (right) stars, both 
showing a low-density, oblate envelope surrounding a 
compact spherical nucleus with a central density almost 
unaltered with respect to that of the parent stars. This 
particular configuration is common to almost all the 
other collisional products as shown in Figure which 
gives column density plots of the ejected stars. Rotation 
as well as asymmetric shape can have a fundamental role 
in the future evolution of the stars and their observable 
characteristics; a star with a rapidly rotating nucleus can 
have its main-sequence life-time considerably extended 
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Fig. 7. — Stars ejected in our simulations after a collision with the companion star. The stars show, typically, an oblate envelope 
surrounding a high density spherical nucleus. Only in simulation C4, where the impact is more "grazing," the collisional product is 
spherically symmetric even in its outermost envelope. 



with respect to their non-rotating counterpart (|Clementl 
1994). Only in run C4 the HVS ejected after the collision 
is slowly spinning and spherical even in its outermost en- 
velope. In this simulation, the stars experience a more 
grazing collision which leads to some envelope-ejection 
but leaves the stars' structure essentially unchanged. 

4.2. Mergers 

Stellar collisions due to either binary evolution or 
dynamical interactions are thought to be the main 
formation channel of blue stragglers in star cluster s 
(|Collier et al.l[l98i iLeonardl fl9M iMateo et al. I fl990h . 
Similar processes have been proposed in the past to ex- 
plain the puzzling presence of the young massive stars 
observed at galactocentric distances of few mpcs, where 
star formation is thought to be strongly inhibited by 
the S MBH tides (jGenzel et all 120031 : lEisenahauer et all 
2005). In paper I, we showed that gravitational encoun- 
ters involving stellar binaries and the SMBH lead, for a 
wide range of orbital parameters, to a stellar collision and 
that among the collisional products, stellar coalescence 
occurs in more than 80% of the cases. In this section we 
study this latter outcome and investigate the properties 
of the resulting stars to clarify whether they would be ex- 
pected to posses features commonly associated with the 
S-star population. 

Table [3] gives the orbital parameters (i.e., eccentricity 
and semimajor-axis) of the merger products in our sim- 
ulations as well as the mass ejected from the binary and 
the fraction of mass captured by the SMBH after the 
first periapsis passage. The table shows that the merger 
remnants lie on a orbit very close to the initial orbit of 
the center of mass of the binary around the black hole, 
implying only a small effect of the mass-loss on the dy- 
namical evolution of the stars. The mass ejected from 
the binary after the first periapsis passage is typically 
larger than that found in collisions that do not end up 
with a merger (see Table [2]) and is of order ~ 10 -2 times 
the initial mass of the binary. In many cases most of the 
mass ejected from the binary during the merger remained 



bound to the SMBH. This is a quite different situation re- 
spect to that found in Section |4~H where approximately 
half of the mass ejected from the binary remained un- 
bound to the black hole; in these previous runs one of 
the two stars is always found on an escaping trajectory 
and, consequently, the debris associated with such a star 
will also tend to escape the SMBH. The last column in 
the table gives the dimensionless spin parameter defined 
in equation ([9]) of the final merger products that show 
very large spins, some of them close to the "break-up" 
value (i.e., J = 1). 

In principle it is possible that, as a consequence of the 
mass loss occurring near periapsis, the resulting merger 
product gains orbital energy and escapes the SMBH (see 
for instance Faber et al. [2005]). Although we do not 
exclude this outcome for a different set of initial condi- 
tions, in our simulations this mechanism does not pro- 
duce HVSs, and in all cases the merger remnant is still on 
a bound orbit around the SMBH. Another more impor- 
tant consideration is that, subsequent to the merger, the 
tidal heating results in some degree of expansion and a 
weakly bound configuration for the merger product. Be- 
cause the tidal radius of the newly formed star is much 
larger than that of its progenitors, the star will succes- 
sively lose more mass with each periapsis passage and 
eventually be torn apart by the SMBH tides. And in 
fact, as it is shown below, this is the final outcome of 
some of our simulations. 

An example of merger is displayed in Figure [5J which 
involves a binary with ao = 0.2 au and equal-mass com- 
ponents of masses M = 3M (run M4). In this run the 
stars collide for the first time after ~ 4 days from the time 
corresponding to the periapsis passage; the internal pe- 
riapsis separation at the first contact is ro/(2R) = 0.86. 
After the first periapsis passage, as consequence of the 
SMBH perturbation, the binary star becomes very ec- 
centric. As the stars move through the circum-binary 
envelope formed during the previous encounters, the or- 
bit gradually circularizes and shrinks. By t « 60 days, 
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TABLE 3 



bame as Table [2] but for stellar mergers after the first periapsis passage. Here the quantities in parentheses refer to the initial orbit of the binary 
center of mass. 



Run 


a 

(au) 


e 


AM h /M h 




AM,/M b 




J 




Ml 


1060(1000) 


0.995(0.995) 


— 2.57 x 10" 


2 


2.57 x 10" 


2 





153 


M2 


984(1000) 


0.999(0.999) 


-1.61 x 10" 


2 


1.26 x 10~ 


2 





167 


M3 


997(1000) 


0.996(0.996) 


-5.07 x 10~ 


2 


2.71 x 10~ 


2 





103 


M4 


1010(1000) 


0.996(0.996) 


-4.26 x 10" 


2 


4.24 x 10" 


2 





259 


M5 


1020(1020) 


0.958(0.963) 


-2.16 x 10" 


2 


1.86 x 10" 


2 





222 


M6 


996(1000) 


0.997(0.997) 


-5.60 x 10~ 


2 


2.93 x 10~ 


2 





0666 


M7 


1180(1000) 


0.997(0.997) 


-2.50 x 10" 


2 


2.37 x 10" 


2 





242 


M8 


1030(1000) 


0.996(0.996) 


-4.00 x 10" 


2 


2.06 x 10" 


2 





249 


M9 


1000(1000) 


0.995(0.995) 


-2.93 x 10~ 


2 


2.65 x 10" 


2 





238 


M10 


1020(1000) 


0.998(0.998) 


-6.09 x 10~ 


2 


5.99 x 10~ 


2 





108 


Mil 


1010(1010) 


0.980(0.982) 


-4.38 x 10~ 


2 


2.51 x 10" 


2 





131 


M12 


1020(1010) 


0.992(0.991) 


-1.80 x 10" 


2 


1.23 x 10" 


2 





181 


M13 


1060(1000) 


0.997(0.994) 


-2.20 x 10" 


2 


1.63 x 10" 


2 
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Fig. 8. — Column density plots for simulation M4 on the X — Z plane. The binary has an internal semimajor-axis ao 
components have masses M = 3Mq. Time t = corresponds to the periapsis passage of the binary external orbit. 
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after approximately 30 collisions, the two stellar nuclei 
merge. The final product has an oblate shape which is a 
common characteristic of all the merger remnants formed 
in our simulations. 

As another example, Figure [9] displays column density 
plots of simulation M13 where a merger occurs between 
a 6 and a 1 M stars. The stars collide after 1.48 days 
from the moment of the closest approach to the SMBH, 
and subsequently merge in the following ~ 1 day. During 
the merger, the high density core of the lower mass star 
rapidly sinks to the center of the companion star. The 
tail-like feature observed at t = 1.86 days in the figure, is 
mostly material coming from the the secondary star that 
loses part of its outermost envelope while sinking to the 



center of the merger remnant. 

Figure [10] shows chemical composition profiles of the 
merger products for runs M4 and M13, after one periap- 
sis passage. In simulation M4, the remnant has a mass 
of roughly 5.9 M , its composition profile is very sim- 
ilar to that of the parent stars (see Figure [2]). Based 
on how long it would take a normal star of that mass 
to evolve to that central hydrogen abundance using the 
TWIN stellar evolution code, we estimate that a "nor- 
mal" 5.9 Mq star would reach a core helium abundance 
Y = 0.35 (assuming Z = 0.02) after - 2Myr. By col- 
liding two 50 Myr old 3M stars, we have effectively 
made a more massive (M ~ 6M ), younger (age ~ 2 
Myr) star. The merger remnant of run Ml 3 shows a pe- 
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Fig. 9. — Column density plots for simulation M13 on the X — Y plane. The binary has an internal semimajor-axis ao — O.lau and its 
components have masses M\ = 6Mq and M2 = IMq. Time t = corresponds to the periapsis passage of the binary external orbit. 



culiar composition profile when compared to a normal 
star, but quite normal for merger products. In Ml 3 (and 
in M12 as well), the low mass star drops to the center 
of the merger product bringing its fresh hydrogen fuel 
along and significantly rejuvenating the core. As a re- 
sult, the maximum He does not occur at the center of 
the merger product. The main reason of the negligible 
amount of hydrodynamic mixing is that the cores of the 
initial stars are very dense and difficult to break even in 
a head-on collision (jLombardi et al. 1I1995L Il996[ ). How- 
ever, we cannot exclude the possibility that other pro- 
cesses occurring on a thermal timescale (as opposed to 
the hydrodynamic timescale) can produce a significant 
degree of mixing in the stars as th ey evolve toward ther- 
mal equilibrium ([Sills et al.lll997| ). 

Another interesting result, plotted in Figure [lOj is 
that a large fraction of the Lithium/Beryllium/Boron 
from the parent stars gets ejected, indicating that a 
significant gap in the abundances of these elements in 
the S-star population might be observational evidence 
of rejuvenation through merger. Reduced atmospheric 
Lithium abundances are for instance observed in blue 
stragglers and can also be a strong indicator of m ixing 
(jHobbs fc Mathieul[l99ll : iPritchet fc Glaspevl[T99Tl ). 

We finally note that as the stars keep orbiting the 
SMBH, their chemical profile and their spinning con- 
figuration will change in time and therefore the states 
displayed in Figures [8j [9] and [10] should be intended as 
not permanent. The evolution will typically lead toward 
smaller spins and a lower Lithium/Beryllium/Boron 
abundances in the stars. We will come back to this point 
below. 

4.3. Clean Ejection of Hypervelocity Stars 



Even when the stars do not collide, tidal torque and 
mass-loss can occur if the periapsis distance of the binary 
center of mass initially lies within the Roche limit of its 
member stars. Similarly to what is done in the previous 
two subsections, we analyze here the first binary- SMBH 
interaction, while we discuss the following evolution of 
the bound stars in the next subsection. Some of the re- 
sults of our SPH simulations, for which there is not a 
direct collision between the stars, are listed in Table [3] 
As expected, for £ > 1 there is no mass-loss, and the 
stars maintain their initial configuration essentially un- 
altered (runs H2 and H3). Conspicuous mass- loss instead 
occurs for runs H6 and H8 in which at least one of the 
stars crosses its tidal radius. Interesting, although the 
usual condition for tidal disruption is well satisfied (i.e., 
A < 1), in both runs the stars are not fully disrupted 
by the SMBH's tidal gravity at the first periapsis pas- 
sage. In the table we also give the final values of the 
spin parameter J that, in general, are found to be a very 
small fraction of the breakup value. Figure [TT1 shows the 
temporal evolution of J and the stellar masses for the 
cases of Table [3] that have the highest value of the final 
spin. The interaction with the SMBH induces a strong 
rotation only in the stars of runs H6 and H8. We con- 
clude that, unless the stars penetrate deeply their tidal 
disruption radius, it seems unlikely that the SMBH tides 
at periapsis alone can produce a significant spin-up of 
the stars. 

As an example, Figure [12] gives column density plots of 
run H8. In this simulation, the primary and secondary 
stars have masses 1 M and 6 M respectively. The pe- 
riapsis of the external orbit (~ 0.8au) is initially inside 
the tidal radius of the 6 M member but it is still outside 
the tidal radius of the 1 M star. At periapsis, the stars 
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Fig. 10. — Upper panels: composition profile of the merger rem- 
nant formed after one periapsis passage in run M4 (see also Figure 
[8]). The merger product of two equal mass stars has a composi- 
tion profile very similar to that of its parent stars. For a ~ 6 M© 
star, a final He abundance in the core of Y = 0.35 corresponds 
to an effective age of ~ 2 Myr. Lower panels: composition profile 
of the merger remnant formed after one periapsis passage in run 
M13 (see also Figure [9]). In this case, the merger product has a 
peculiar profile if compared to a "normal" star. Its core is strongly 
hydrogen-enriched as a consequence of the low-He fluid transported 
by the low mass star along with it to the center. 



are squeezed by the SMBH's tidal gravity. In the process 
the primary star losses a large fraction of its initial mass 
(~ 3.4 M ), while the secondary loses only ~ 0.03 M . 
After the interaction with the SMBH the binary is bro- 
ken apart and the lightest member becomes an HVS. Be- 
cause of tidal heating during the periapsis passage, the 
stars are perturbed from their thermal equilibrium state 
and their radii are somewhat enlarged with respect to a 




Fig. 11. — Like Fig. [5] but for some of the simulations in which 
there is neither a stellar collision nor a merger. 



normal main-sequence star of the same mass. 

4.4. The Bound Population 

After the initial encounter between the binary and the 
SMBH, one star remains in a bound orbit around the 
SMBH in all of our simulations. Like the orbit of the ini- 
tial binary about the SMBH, the orbit of such a bound 
star is highly eccentric: 0.96 < e < 1. In those cases 
in which the bound star is a merger product (runs Ml 
through M13), the semimajor-axis a very nearly equals 
the semimajor-axis of the initial binary about the SMBH: 
a ~ lOOOau, corresponding to an orbital period of about 
16 years. We note that these orbital periods are compa- 
rable to those of the S-Stars in the GC. In those cases 
in which a HVS star is ejected (runs CI through C5 and 
runs HI through H8), the ejection energy comes at the 
expense of the orbital energy of the bound star, which 
consequently has a somewhat smaller semimajor-axis: 
lOOau < a < 700au, corresponding to orbital periods 
of 0.5 to 9 years. 

4.4.1. Tidal Stripping 

The most significant hydrodynamic effects occur near 
the periapsis, where induced collisions and mergers are 
most likely to occur and where tidal stripping is at its 
greatest. 

We find that the periapsis separation of a bound star 
remains remarkably constant from one orbit to the next, 
even when there is significant mass loss due to Roche 
lobe overflow at periapse. As an example, consider the 
run CI in which the initial stars have the dimensionless 
Roche lobe parameter £ = 0.67. As expected for ( < 1, 
the bound star does indeed lose mass each time it sweeps 
past the black hole. The gradual decrease of the mass 
M of the bound star can be seen in the top panels of 
Figure [131 The mass 5M lost per orbit, shown by the 
star symbols in the top panel, increases with each or- 
bit, until after 96 periapsis passages the star has been 
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Fig. 12. — Column density plots for simulation H8 on the X — Z plane. The binary has an internal semimajor-axis ao = 0.1 au and its 
components have masses Mi = 6Mq and M2 = IMq. Time t = corresponds to the periapsis passage of the binary external orbit. The 
first four panels are centered to center of mass of the binary, while in the two bottom-right panels the origin is the center of mass of either 
the captured (left) or ejected star (right). 



completely disrupted. We note from the bottom panel 
that the periapsis separation r per = a(l — e) is nearly 
unchanged during this entire process: in this and other 
cases, we find the bound star returns to the nearly same 
relative separation from the black hole regardless of mass 
loss. The apoapsis separation d = a(l + e) is also some- 
what constant, although it decreases at late times when 
the mass loss is greatest and strong tidal effects remove 
energy from the orbit. 

As another example, the lower panels of Figure [13] give 
the evolution of the merger remnant formed in run M5. 
In this case most of the mass-loss occurs during the first 
periapsis passages where very high entropy material is 
removed from the outer layers of the star that responds 
by reducing its radius. In the following evolution, mass- 
loss essentially ceases. We note that merger products 
have a very non-uniform density profile characterized by 
a extended low density envelope and a dense central re- 
gion. Subsequent passages of the star by the SMBH will 
therefore cause the depletion of the outermost stellar re- 
gion, unveiling its hot central core. An example of this 
phenomenon is shown in Figure [14J where we plot col- 
umn density plots for the merger remnant of run M7. 
After about twenty orbits mass loss stops and the enve- 
lope has been completely removed. Similar mechanisms, 
involving tidal stripping suffered by late-type giants dur- 
ing close passages around a intermediate massive b lack 
hole, have been invoked in the past (jMiocchi II2QQ7I ) to 
explain the extreme ho rizontal-branch st ars observed in 
some globular clusters (jRich et al. Ill997f ). 

We find that the dimensionless parameter ( is strongly 
correlated with whether and how quickly a bound star 



loses mass through Roche lobe overflow. For example in 
run H2, the bound star has ( > 1 and does not experience 
a collision or merger that would change £: it consequently 
continues to orbit the SMBH without ever suffering any 
mass loss. In all the cases with £ < 1, the bound star 
is ultimately destroyed after repeated episodes of Roche 
lobe overflow, with smaller values of ( generally corre- 
sponding to fewer orbits before disruption. Figure [T5l 
shows that the number N p of periapsis passages before 
disruption grows exponentially with the initial (. In ad- 
dition, this number of passages depends only weakly on 
whether the interaction type is a collision (red crosses), 
merger (black triangles), or clean ejection of a HVS (blue 
circles). 

The rightmost data point in Figure [15j corresponding 
to run C2, deserves some discussion. In this case, £i = 
(2 = 1.19 > 1, so that neither binary component would 
lose mass if it were not for the collision induced on the 
first periapsis passage. This collision both increases the 
radius and slightly decreases the mass of the bound star, 
effectively decreasing its £ parameter to a value below 1. 
Thus, on the subsequent passage past the black hole, the 
star loses more mass, now due to Roche lobe overflow. 
The response of this particular star to mass loss is that its 
radius remains roughly constant. From equation the 
( parameter then stays below 1 and slowly decreases as 
the mass ratio q decreases with each successive passage. 
Ultimately, after nearly 600 periapsis passages, the star 
is completely pulled apart. 

Similarly to the stars from run C2, the 6 M star in run 
M13 has £ = 1.22 > 1. This star indeed makes the first 
periapsis passage without immediately losing any mass; 
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Fig. 13. — Evolution versus time (and number of orbits) in the 
runs CI and M5. From the top to the bottom panel: mass 8M m 
gained per orbit by the SMBH (filled circles) and mass SM lost per 
orbit by the stars (star symbols), cumulative mass AM. bound to 
the SMBH, mass M of the bound star, apoapsis d of the bound 
star, periapsis r per of the bound star. 



however, while the binary recedes away from the SMBH, 
the merger causes mass ejection. The resulting 6.84 M 
merger product is large enough that ( drops below 1, 
and, on the subsequent periapsis passages, mass is lost 
through Roche lobe overflow. As a result of shedding 
its high entropy outer layers, the merger product shrinks 
sufficiently that ( is pushed back toward ( « 1. By 
comparing runs C2 and M13, we conclude that the fate 
of binary stars with £ > 1 depends not simply on the 
initial values of £ but also on the type of their interaction 
and the response of the bound star to mass loss. 



In several of our simulations, merger products formed 
from stars with £ > 1 are large enough that £ drops be- 
low 1 and at least some mass is lost due to Roche lobe 
overflow on the second and later periapsis passages (runs 
Ml, M3, M4, M5, M7, M8, M9, Mil, M12, and M13). 
Because shock heating is preferent ially distributed to th e 
outer layers of a merger product (jLombardi et al. 2002), 
this Roche lobe overflow always strips away very high 
entropy material and the product responds by decreas- 
ing its radius. In this way, the ( parameter gradually 
approaches a value « 1, corresponding to an eccentric 
semidetached binary consisting of the bound star and 
the SMBH. In our SPH simulations of such cases, we 
typically follow the dynamics for several hundred orbits, 
without seeing an appreciable decrease in the mass of the 
bound star: indeed at late times the mass loss typically 
fluctuates between and 2 SPH particles per periapsis 
passage. Such situations necessarily challenge the mass 
resolution limit of our simulations, and it is difficult to 
say whether such small levels of mass loss are physically 
meaningful or simply a numerical artifact. In any case, 
in nature, thermal relaxation in the outermost layers of 
such a merger product would tend to retract it inside of 
its Roche lobe and stabilize the star against further mass 
loss. 

To better understand the effects of the numerical res- 
olution, we vary the number of particles used to model 
several of the scenarios. The results for scenario M7 are 
summarized in Table [5j where we list the mass, eccentric- 
ity, and semimajor axis of the bound star after 25 orbits. 
We find a good agreement of results at all resolutions 
tested and a convergence of these results as the number 
of particles N is increased up to the value used in this 
paper 4 x 10 4 ). In particular, the final orbital data 
have converged to within ~ 0.02%. 

We also extend our resolution study to cases in which 
the bound star is ultimately disrupted. We find that, 
for various particle numbers from N ~ 5 x 10 3 up to 
8 x 10 4 , the simulations of the same initial conditions all 
behave very similarly for at least the first ~ 100 orbits 
around the SMBH. The top three frames of Figure [16j 
for example, demonstrate this consistency for the mass 
M, eccentricity e, and semimajor axis a of the bound 
star after 50 orbits. For situations in which the star 
orbits the SMBH more than ~ 100 times, the simulations 
at various resolutions diverge at late times, with higher 
resolutions simulations in which there is a collision or 
merger requiring more periapsis passages to disrupt the 
bound star (see C2 and M6 the bottom frame of Fig.[T6|). 

In addition to the scenarios shown in Figure [16j we also 
study resolution effects in simulations of several other 
cases in which the bound star is ultimately disrupted 
(specifically C3, HI, H3, H4, H6, and M2), with the par- 
ticle number N varying from 5 x 10 3 up to 4 x 10 4 . 

As with the Figure [16] data, the consistency of the re- 
sults is again very good for N p < 100. For example, 
for case C4, each of four such simulations predict that it 
would take somewhere in the range of 89 to 91 periastron 
passages to completely disrupt the bound star. Further- 
more, for case C3 all simulations predict that it takes 8 
periastron passages to disrupt the bound star, in case H4 
all simulations predict that it takes 10 passages, and in 
case H6 all simulations predict that it takes 2 passages. 




Fig. 14. — Evolution of the merger remnant formed in run M7 after each periapsis passage when the star sets to hydrostatic equilibrium 
and until mass loss ceases. Time increases from left to right and form top to bottom. Initially the merger remnant has a large low-density 
envelope that is completely removed after several obits around the SMBH. 



We conclude that the particle number employed in this 
paper is sufficient to model accurately the evolution for 
at least - 100 orbits around the SMBH. 

Finally, as an illustrative example, Figure [T71 gives col- 
umn density plots for run M6 after 50 periapsis passages 
and for simulations with different number of particles. 

4.4.2. Internal Structure 

Figure [18] shows the composition profiles as a function 
of enclosed mass fraction m/M for the bound star in 
cases M4 and M13. Here, m is the mass enclosed within 
an isodensity surface and M is the total bound mass. The 
dotted curves show the profiles after two periapse pas- 
sages, while the solid curves show the same profiles once 
the bound star has effectively reached a steady state. 
Mass loss experienced during multiple passages removes 
the outer layers of the star, decreasing the bound mass 
M and causing the composition profiles to shift slightly 
to larger enclosed mass fractions m/M. 

We note that the helium profile for Ml 3 is qualita- 
tively sim il ar to that of the case G merger product in 
ISills et all (fl997h ( see Fig. 2 in that paper): both have a 
maximum helium abundance at an intermediate radius 
inside the star. In both case G and our M13, the strange 
helium profile is caused by a low mass star sinking to 
the center of the collision product and displacing the he- 
lium rich fluid outward. The stellar tr ack for the case G 
product is shown in Figures 4 and 6 of Sill s et al.l (fl997h . 



In Figure 6, we see that, on the main sequence, the case 
G product is somewhat bluer and brighter than a nor- 
mal main sequence star of the same mass. In Figure 
4, we see that the case G product is a little bluer and 
brighter than a different collision product (case J) with 
basically the same mass but without the dense hydrogen 
core. It is the increased helium content in the stellar in- 
terior that makes the opacity lower (compared to other 
main sequence stars of the same mass) and thus bluer 
and brighter. So, by analogy, our M13 merger product 
would be a little bluer and brighter than a normal main- 
sequence star of the same mass. 

The elements Li, Be, and B are potentially interesting 
observational indicators of the history of a star. These el- 
ements burn at temperatures of about 2.5 x 10 6 , 3.5 x 10 6 , 
and 5 x 10 6 K, respectively, and therefore can exist only 
in thin outer layers of the parent stars. During a dynam- 
ical interaction, these elements can be removed either by 
ejecting the outer layers of the star or by redistributing 
them to an environment too hot for their long term ex- 
istence. Although Be, B, and Li can exist after the first 
periapsis passage (see Fig. [TP]), the effect of multiple pas- 
sages is typically to remove these elements completely. 
Although there are cases where B still exists in the final 
bound star, it is always severely depleted. For example, 
in run Ml 3 the B level at the surface of the final product 
is only ~3% of the surface value in the 6 M parent star 
from which it originated. 
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Fig. 15. — The number of periapsis passages past the black hole 
needed to disrupt the bound star versus the initial dimensionless 
Roche lobe parameter £i. The different data points represent the 
scenarios in which the bound star suffers a collision (red crosses: 
runs CI, C2, C3, and C4), is formed in a merger (black triangles: 
runs M2, M6, and M10), or is cleanly separated from the HVS (blue 
circles: runs HI, H3, H4, H5, H6, H7, and H8). Those collision 
and merger data with N p > 100 likely underestimate N p due to 
resolution effects and are best considered as lower limits (see text). 
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Fig. 16. — The mass M in solar masses, eccentricity e, and semi- 
major axis a of the bound star, all after 50 periapsis passages, 
versus total particle number N for several representative scenar- 
ios in which the bound star is ultimately disrupted: CI (red), C2 
(green), C4 (blue), H5 (cyan), M6 (magenta), and M10 (black). 
Also shown, in the bottom frame, is the number N p of periapsis 
passages needed to completely disrupt the bound star. For a given 
scenario, note the consistency of the data for M, a, and e after 50 
orbits, as well as for N v in cases with N v < 100. 



In Table El we summarize some properties of the bound 
stars that survive in our SPH simulations (i.e., that are 
not ultimately disrupted by the SMBH). The "number 
of passages" represents the number of periapse passages 
before the mass loss effectively shuts off, which we de- 
fined as having 2 or fewer SPH particles ejected. The 
central hydrogen abundance is given by X c , which al- 
ways equals the central hydrogen abundance of the low- 
est mass parent star. We also list the effective age of 
the bound star based on its mass and central hydrogen 
abundance. We evaluate this effective age, based on how 
long it would take a normal star of that mass to evolve to 
that central hydrogen abundance using the TWIN stel- 
lar evolution code. It is known that the contraction of a 
merger product to the main-sequence is very similar to 
the contraction of a pre-main-sequence star to the main- 
sequence. In the latter case, the most important variable 
is the mass (for a given hydrogen abundance). In the 
former case, the two important variab les are essentially 
mass and central hydrogen abundance ( Sill s fc Lombardil 
fl997h . In runs C5 and H2 the bound star is only a slightly 
perturbed version of one of the binary components, and 
the ages in these cases is the same 50 Myr as that com- 
ponent. For mergers of two 3M stars, the effective age 
of the merger product is in the range of 14 to 22 Myr. 
For mergers of two 6M stars, the effective age is in 
the range of 6 to 9 Myr. Mergers of unequal mass stars 
(M12 and M13) also significantly rejuvenate a star: for 
example in M13, the sinking of the 1 solar mass star to 
the center of the merger product essentially resets the 
nuclear clock to only 0.3Myr after the ZAMS. 

In the last three columns of Table [6] we list central 
temperature T c , internal energy U and thermal timescale 
^thermal of surviving bound stars. The central tempera- 
ture is defined as the temperature in the star where the 
density is highest and therefore is not always the tem- 
perature of the highest temperature SPH particle. The 
central temperature T c of all the stars is large enough to 
sustain nuclear burning in the core. The global thermal 
time scale can be estimated as 



^thermal 



U/(L), 



(10) 



where (L) is a mass weighted average of the luminosity 
L throughout the entire star. 

To calculate L, we take advantage of the fact that 
the parent stars are massive enough to be fully radia- 
tive. In addition, shock heating prevents any convective 
zones from existing in a newly formed merger product. 
Thus, we obtain the luminosity L exiting a closed sur- 
face by the integral L = § F • da, where da is an area 
element on the surface and the diffusive radiative flux 
F = — 4acT 3 VT/(3ft/>). Here a is the radiation constant, 
c is the speed of light, and n is the opacity. The sur- 
face integral is easily converted to a volume integral by 
the divergence theorem. The result, L = J V • FdV, is 
straightforward to estimate in SPH: 



Vfli 



(11) 



where the sum is over only those particles positioned in- 
side the surface under consideration. Because SPH calcu- 
lations cannot properly resolve the photosphere, equation 
(pTj) cannot be used to give a reliable total luminosity. 
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Fig. 17. — Column density plots for run M6 after 50 periapsis passages and using different total number of particles (5k, 10k, 80k). 



However, equation (fTT|) does allow us to study the lumi- 
nosity profile throughout the bulk of the system. 

As an example, in Figure [T9j we show a more detailed 
look at the interior of the final bound stars in runs M4 
and C5. From top to bottom, we give the luminosity L, 
temperature T, and radius r as a function of enclosed 
mass m. To evaluate the luminosity profile, we use equa- 
tion (JTTJ on each SPH particle, summing over particles 
of larger density and calculating the opacity k from the 
OPAL tables. Our merger and collision products typi- 
cally achieve a maximum luminosity in their outermost 
layers that is comparable to the Eddington luminosity 

L e dd = — M~3.8xlO 4 L -^- , (12) 

although such a high luminosity would diminish rapidly 
as the star contracts to the main-sequence in a time 
thermal usually of order ~ O.lMyr. 

5. DISCUSSION 

The observed rotation rates of HVSs may give i mpor - 
tant clues to their formational history. iHansenl (j2007l ) 
has proposed that, as a consequence of tidal locking 
in close binaries, HVSs ejected by the Hills mechanism 
should rotate system atically slower than field s t ars o f 
similar spectral type. lLopez-Morales fc Bonanosl ([2008) 
found that the late B-type star HVS 8 has a rota- 
tional velocity of — 260 kms -1 , more typical of sin- 
gle B-type stars and therefore seemingly contrary to 
the hypothesis of a binary origin for this star. In or- 
der to explain the observations, other ejection mecha- 
nisms have been invoked, such as ejection by a close 
encounter with a massive black hole binary or with a 
stellar black hole orb i ting the ga l actic center S MBH 
(lYu fc Tremaind 120031: [Le vin 20061: iSesana et all 120061 : 
lLockmann. fc Bau mgardt 2QQ8h . However, it has been 
note that a larger statistic would be certainly required in 
order for the rotation to be us ed as a signatu re for the 
origin of HVSs and/or S-stars (jPerets Il2009af ). Further- 
more, in this paper we have shown that there are two 
other potentially important ways with which the stars 
can somewhat increase their rotation even in the binary 
disruption scenario: tidal torque by the SMBH at periap- 
sis (if the stars enter within their tidal disruption radius) 
and/or a collision between the two binary members. 

With the help of simplifying approximations, we are 



able to relate the rotational parameter J calculated in 
our simulations to the observable rotational velocity v 
after that star has thermally relaxed back to the main 
sequence. In particular, we approximate that the star 
rotates rigidly, that its rotation does not drastically af- 
fect its structure, and that the rotational parameter J is 
conserved during relaxation. Using L = Iv/R = c\MRv 
and E = —C2GM 2 /R in equation ([9]), we obtain J = 

c 1 c\ /2 v{R/{GM)) 1 ^. Clearly c x and C2 are simply nu- 
merical coefficients related to the moment of inertia / 
and total energy E of a star, respectively. For B-type 

1/2 

main sequence stars, we find c\c 2 ~ 0.04 to 0.05 and 
R/M -0.5 to 0.8 R /M using models from the TWIN 
stellar evolution code. Solving for v in terms of J, we 
find 

v ~ 1.2 x 10 4 Jkms _1 , (13) 

accurate to within — 30% for most B-type main sequence 
stars. Given the J values of ejected stars in our simula- 
tions (see Fig. [5] and Table HJ, we estimate from equation 
([T3|) that the post-relaxation rotational velocity v can 
be as large as —400 or 500 km s _1 for HVS stars (con- 
sider runs C3 and H6). We therefore conclude that the 
rotation of the star HVS 8, for example, is completely 
consistent with a binary origin. 

The J values for unmerged stars in our simulations that 
are bound to the SMBH after the first periapsis passage 
indicate, via equation (|T3|) . that the post-relaxation ro- 
tational velocity v would be typically < 400km s _1 but 
could be as large as — 1000km s _1 . These large rotation 
velocities, however, correspond to stars that penetrate 
deeply within their tidal disruption radius (e.g. runs H6 
and H8) and therefore are eventually destroyed after sev- 
eral orbits. Merger products obtain even larger spins af- 
ter the first periapsis passage (see Table [3|). However, 
as the merger product keeps orbiting the SMBH, J de- 
creases as mass gets pulled off the outside of the star, 
where the specific angular momentum is greatest. The 
bound stars that survive to orbit the SMBH end our sim- 
ulations with J < 0.017 (see Tabled]), corresponding to 
post-relaxation rotational velocities v < 200km s _1 . As 
our simulations began with irrotational stars, the actual 
final rotational velocity of a (bound) star could be larger 
or smaller if the parent stars had significant spin, de- 
pending on the orientation of the spin axis with respect 
to the external orbital plane and/or the plane on which 
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TABLE 4 



Same as table [2] but for simulations in which the stars do not collide and one member becomes a HVS. Jcaptured gives the spin of the stars that 
remain bound to the SMBH, while J e jected refers to the ejected stars. 



Run 


(km/s) 


a 

(au) 


e 


AM h /M h 




AM. /M b 




•^captured (^ejected ) 




HI 


3883(3910) 


160(160) 


0.991(0.991) 


-2.230 x 


io- 4 


1.221 x 10" 


-4 


0.211 x 10" 


• 2 (0.188 x 10" 


- 2 ) 


H2 


2026(2160) 


309(308) 


0.983(0.974) 












0.268 x 10" 


' 3 (0.238 x 10" 


- 3 ) 


H3 


1082(1100) 


423(421) 


0.991(0.983) 












0.153 x 10" 


3 (0.729 x 10" 


- 4 ) 


H4 


4221(4252) 


142(142) 


0.995(0.990) 


-1.019 x 


10~ 2 


5.401 x 10" 


-3 


0.194 x 10" 


1 (0.147 x 10" 


- 1 ) 


H5 


2121(2155) 


306(307) 


0.974(0.993) 


-8.454 x 


10~ 5 


4.759 x 10" 


-5 


0.149 x 10" 


' 2 (0.168 x 10" 


- 2 ) 


H6 


1637(947.6) 


174(444) 


0.988(0.999) 


-0.9914 




0.9914 




0.445 x 10" 


1 (0.379 x 10" 


- 1 ) 


H7 


2570(2574) 


677(674) 


0.983(0.978) 


-1.335 x 10~ 2 


6.911 x 10" 


-3 


0.223 x 10" 


1 (0.348 x 10" 


- 3 ) 


H8 


2494(2497) 


619(619) 


0.991(0.990) 


-0.4835 




0.2458 




0.845 x 10" 


1 (0.274 x 10" 


- 1 ) 



TABLE 5 



Resolution study for scenario M7. The particle number is given by N, 
while the mass, orbital eccentricity and semimajor-axis of the bound 
star after 25 orbits around the SMBH are given by M, e, and a respec- 
tively. 



N 


M 


e 


a 




(M ) 




(au) 


4,957 


0.9972 


1147.1 


8.468 


9,889 


0.9973 


1146.7 


8.410 


19,933 


0.9974 


1146.2 


8.392 


39,877 


0.9973 


1146.2 


8.390 



the collision occurs. If the spin axis is aligned with the 
angular momentum of the external orbit, the initial spin 
will sum up with that acquired due to the tidal torque 
from the SMBH. A larger spin will also result from a 
collision, if the angular momentum of the inner binary 
is initially aligned with the spin axis of the stars. We 
stress here that, in general, the effect of an initial spin 
on the final rotation of the stars can be complicated, for 
this reason we decided to ignore initial rotation and be- 
gin with binary components irrotational in the inertial 
frame, which allows us to more easily measure any rota- 
tion imparted during the subsequent interaction. 

Deep near-IR observations of the GC show that the S- 
stars are B0-B9 main-sequence stars with rotational ve- 
lociti es similar to thos e of field stars of the same spectral 
type (| Alexanderl [2005) . Our final bound stars therefore 
have properties very similar to those of the S-stars: their 
masses qualify them as spectral type B main-sequence 
stars, and their post-relaxation rotational speeds are of 
the correct general magnitude. For example, the rota- 
tional speeds of our fastest rotators are con sistent with 
the 2 20 ±40 km s _1 value for the S-star SO-2 (|Ghez et al.l 
2003). However, we note that the orbital eccentricities 
of our bound stars (0.96 < e < 1) are larger than that 
of SO-2 (e ~ 0.87 ) , as similarl y found in simulations by 
IGinsburg fc Loebl (|2006D and [Hansen) (|2007D . 

For tidal torque from the black hole to have a signifi- 
cant effect on stellar rotation, the stars should enter deep 
into their disruption zone (i.e., r per < ft). 

When no collision occurs between the components of 
a binary, the distance of closest approach of the two 
stars to the SMBH typically changes little due to the 
encounter. In such cases, a necessary condition for sig- 
nificant spin-up is that the binary itself be on an orbit 



that passes within ~ r t of the SMBH. This implies in 
turn that the fractional change in orbital angular mo- 
mentum with respect to the SMBH, per orbit, be of or- 
der unity. This condition is satisfied in the so-called "full 
loss cone" regime, which, in a galaxy like the Milky Way, 
extends inward to ~ 0.2 times the SMBH influence ra- 
dius, or to r « 0.5pc (e.g. Wang & Merritt 2004). Q 
Inside this region, which is the region of interest for the 
current study, evolution onto loss-cone orbits is diffu- 
sive, and most binaries would be tidally disrupted before 
finding themselves on orbits that intersect ~ r t . We 
note however that in the "massive-perturber scenario" 
the apoapsis distance of the binary is of the order of a 
parsec (i.e., > 0.5pc) and therefore the fractional change 
in orbital angular momentum with respect to the SMBH, 
per orbit, can be of sufficient to put the binaries on a 
trajectory that passes within r t of the SMBH. Even in- 
side 0.5p c, other dynamical proce sses like resonant re- 
laxation (Ra uch fc Tremaind 119961), scattering from an 
intermediate mass black hole (jMerritt et al. I l2009f) or 
perha ps eccentric instability in a disc (|Madigan et al. I 
l2008h can produce larger changes in orbital angular mo- 
mentum than in the case of two-body relaxation alone. 
On the other hand tidal spin-up is not expected to be 
very efficient because it is important only for the narrow 
range of periapses: \r t < r < r t (for r < |r*£ the star 
is fully disrupted; for r > r t tidal torque is small). As 
a consequence of the previous condition, the ejected star 
will lose a large fraction of its mass. An observational 
indicator of the history of the star would be, even in this 
case, a deficit in the abundances of light elements (such 
as Lithium) that can exist only in thin outer layers of the 
parent star and that are typically removed by the tidal 
interaction with the SMBH. 

One of the main arguments against rejuvenation of 
the S-stars throu gh merger i s the apparent "normality" 
of their spectra (jFigerl |2008[ ) . We note, however, that 
the envelope of our merger products will not look sig- 
nificantly different than that of normal stars (compare 
the right edge of the plots in Figure 15 with the right 
edge of Figure 2). In fact, if the parent stars are of equal 
mass, the merger product will have very normal profiles 
throughout the star. If the parent stars are of signifi- 
cantly different mass, then the profiles are more peculiar 
and one should worry about how this affects the stellar 
evolution. The main effect would probably be to change 

1 Assuming ap~r~ 2 density cusp. 
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TABLE 6 



Some properties of the bound stars that survive in our SPH simulations. "Number of Passages" gives the number of SMBH-star encounters 
before mass-loss ceases. Mass, orbital eccentricity and semimajor-axis of the star are given by M, e and a respectively. The value of J is the 
final dimensionless spin parameter, while X c is the central hydrogen abundance. The corresponding effective age is also listed. In the last three 
columns, we give central temperature T c (the temperature where the density is highest), internal energy U, and thermal timescale ithermai- 



Run 


Number of Passages 


M 


e 


a 


J 


x c 


Effective Age 


T c 


U 


^ther 






(M ) 




(au) 






(Myr) 


xlO 6 (K) 


XlO 50 (erg) 


(My 


C5 


8 


2.998 


0.974 


212 


0.009 


0.63 


50 


22 


0.141 


0.4 


H2 


2 


3.000 


0.978 


316 


0.000 


0.63 


50 


23 


0.148 


1 


Ml 


30 


4.83 


0.995 


1010 


0.009 


0.63 


16 


31 


0.326 


0.1 


M3 


30 


4.94 


0.996 


996 


0.010 


0.63 


16 


31 


0.321 


0.2 


M4 


18 


4.27 


0.996 


1000 


0.009 


0.63 


22 


32 


0.292 


0.2 


M5 


22 


5.15 


0.962 


1020 


0.014 


0.63 


14 


32 


0.309 


0.05 


M7 


25 


8.39 


0.997 


973 


0.016 


0.56 


9 


39 


0.710 


0.1 


M8 


90 


10.7 


0.996 


1030 


0.014 


0.56 


6 


39 


0.944 


0.07 


M9 


27 


8.74 


0.995 


1010 


0.010 


0.56 


9 


39 


0.736 


0.06 


Mil 


14 


10.1 


0.980 


1010 


0.005 


0.56 


7 


36 


0.793 


0.03 


M12 


27 


7.27 


0.992 


1020 


0.009 


0.63 


7 


32 


0.269 


0.03 


M13 


79 


6.32 


0.997 


1060 


0.017 


0.70 


0.3 


18 


0.391 


0.1 



the opacity and therefore shift slightly the color and lu- 
minosity. But, unless significant mixing is induced, the 
chemical composition of the outermost envelope will re- 
main similar to that of the higher mass star (with the 
possible exception of Li, Be, and B levels: see §4.2|) . 

The tidal disruption of a star passing close enough 
by a SMBH to enter its tidal radius, produces a lumi- 
nous UV/X-ray flare of radiation as the bound s tellar 
gas f alls back onto the black hole and is accreted (Rees 
1988). Tidal flares are of great interest because they 
can probe the presence of SMBHs in galaxies with oth- 
erwise no evidence of an active nucleus and can be used 
to measure the mass and spin of the cen t ral bl ack hole 
(jKomossa et all l2QQ4t iGezari et all I2QQ8L l2QQ9h . Com- 
putations of the tidal disruption of stars have been per- 
formed by several studies in the past, with the aim of un- 
derstan d ing th e observationa l signa t ures of these events 
(lUlmerl [T^g; [Bo^anmdll |2QQ5 iGomboc Cadez I 
2005; Lodato et al. 2009; Strabbe & Quataert 2009; 
Guillochon et al.l 120091 : [Kasen & Ramirez-Ruiz 2010f ). 
We note that there are many important scenarios in our 
paper that have been so far almost completely ignored: 
(i) multiple passages by the same star; (ii) merger of two 
stars resulting in disruption due to the increased size; 
(hi) partial tidal disruptions. 

Predicting the radiative effects of multiple passages of 
a star by a SMBH is outside the scope of the present 
work. But, it seems likely that the light curve result- 
ing from these repeated tidal events might show a se- 
ries of small peaks, separated roughly by the orbital pe- 
riod, before finally producing the large peak that is ob- 
served as the "tidal disruption." Assuming the star on a 
parabolic trajectory, after the star-SMBH encounter, the 
most bound material moves on a orbit with semi-major 
axis a = ^r^ er /R and returns to periapsis after a time 



(GM.) 1/2 (2Rf/ 2 
that it is also the time when the flare starts, while the 



peak return rate occ urs at t ~ 1.5to (|Evans fc Kochanekl 
fl989l : lLi et aHl2002l ). 

As previously discussed, merger products are very 
large, so it is easy to strip off lots of mass during the early 
periapsis passages, while at later time the mass loss of- 
ten ceases. The light curve resulting from these repeated 
tidal events will eventually show a series of small peaks 
of declining intensity (see Figure 12). The result of the 
repeated (partial-)tidal disruption of a star with a large 
envelope (e.g., late- type giants) will show a similar light 
curve. In collisions without mergers, there is some ex- 
pansion in size but it is not as dramatic as in a merger. 
So it is not until late times that there is significant mass 
overflowing the Roche lobe. The light curve will show 
peaks of increasing intensity until the last brightest flare 
produced by the full tidal disruption of the star. 

In future work, one could model the material that be- 
comes bound to the black hole more carefully. In par- 
ticular, it would be interesting to identify possible sig- 
natures of interactions between the bound star and the 
accretion torus left behind from previous periapse pas- 
sages. Quasi-periodic emission may be detected if X-ray 
flares arise ever y time the star crosses the torus plane 
dDai et al.ll20T(l . 

6. SUMMARY 

In this paper, we carried out hydrodynamic simulations 
of binary stars in orbit about the supermassive black hole 
(SMBH) at the Galactic center. In the TV-body simula- 
tions of Paper I, we assigned physical sizes to stars based 
on a simple mass-radius relation and predicted which bi- 
naries would merge, i.e., undergo a collision with relative 
velocity less than escape velocity. The fluid simulations 
presented in this paper were found to be quite consistent 
with the TV-body simulations, in the sense that when the 
latter predicted a stellar merger, the fluid stars typically 
merged as well. The merger rates presented in that pa- 
per are therefore confirmed by the present work. The 
principal, new findings of our work are summarized here. 

1 The central temperature of the merged stars in all 
our simulations is large enough that there would 
still be nuclear burning in the stellar core. How- 
ever, mergers tend to "rejuvenate" stars, in the 
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Fig. 18. — Upper panels: composition profile of the merger 
product in run M4 after two (dotted curves) and eighteen (solid 
curves) periapsis passages, corresponding to masses M = 4.75 Mq 
and M = 4.27 Mq, respectively. Lower panels: composition pro- 
file of the merger product in run M13 after two (dotted curve) 
and 76 (solid curves) periapsis passages, corresponding to masses 
M = 6.46 M and M = 6.32 M , respectively. 



sense that the lower mass star sinks to the cen- 
ter of the merger remnant, effectively resetting the 
nuclear clock of the merger product to the zero-age 
main sequence. Even the products of equal-mass 
mergers exhibit significant rejunvenation. 

2 Mass loss during collisions is generally small, but 
large fractional mass loss can occur when one or 
both stars is tidally perturbed by the SMBH; if the 
two stars merge, a temporarily more extended ob- 
ject is formed, further enhancing the mass loss rate. 




Fig. 19. — From top to bottom: luminosity L, temperature T, 
and radius r as a function of enclosed mass m for the final bound 
stars of runs M4 (upper panels) and C5 (lower panels). L, r and 
m are in solar units, while T is in Kelvin. 



When the merger product has a distance of clos- 
est approach to the SMBH smaller than its Roche 
limit, total disruption always occurs after repeated 
periapse passages. Repeated tidal flares, separated 
by roughly the orbital period, are predicted to pre- 
cede the disruption. 

3 In stellar mergers, elements that can exist only in 
the outermost envelope of the stars such as Li, Be 
and B are severely depleted, due primarily to tidal 
truncation by the SMBH in subsequent periapse 
passages. However, the envelopes of the merger 
products do not otherwise differ significantly from 
those of the parent stars. 
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TABLE 7 

HVSs (v e j > 1000 kms -1 ), Collision and Merger frequency(%) 



M 2 (M ) 


d(pc) 


HVSs (total) 


HVSs (primary) 


HVSs (secondary) 


Collisions 


Mergers 


1 


0.01 


33. 6 


1.29 


32.3 


7.74 


7.10 


3 


0.01 


52.9 


25.5 


27.4 


9.35 


8.71 


1 


0.1 


62.5 


31.3 


31.2 


6.77 


6.45 


3 


0.1 


59.3 


29.3 


30.0 


11.6 


11.3 



We finally stress that SPH calculations neglect radia- 
tive and heat transport, and therefore can follow the sys- 
tem only over hydrodynamical timescales that are typ- 
ically of the order of a few hours. For these reasons, 
in this paper, we were able to discuss the relaxed struc- 
ture of merger products only qualitatively and the re- 
laxation time only in order of magnitude. In a subse- 
quent paper, we plan to present the results of stellar 
evolution calculations that can follow the evolution of 
the SPH merger products over much longer, thermal and 
nuclear timescales and determine their track in a color- 
magnitude diagram. These calculations will allow us to 
compare the observable properties of our models with the 
properties of stars observed at the galactic center. 
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APPENDIX 

UNEQUAL MASS BINARIES: TV-BODY SIMULATIONS 

In the following, we briefly present the results of new TV-body simulations of unequal-m ass binaries passing by Sgr 
A*. The integrations are performed by using the high accuracy TV-body code ARCHAIN (Mik kola fc Merritt 1 12008. 
2006). The binaries are initially placed at a distance d = 0.1 — O.Olpc from the SMBH with a purely tangential velocity 
corresponding to periapses between the tidal disruption radius of the secondary star and r D t- The primary star has 
mass Mi = 6M , while the mass of the secondary is either M2 = 1M or 3M . We adopt ao = O.lau for the internal 
semimajor-axis of the binaries. In all the simulations the final integration time is fixed to one orbital period of the 
binary orbit around the SMBH. In total we perform 1200 simulations. 

As expected, we find a systematically larger ejection velocity for the less massive star (|Yu fc Tre maine 2003). For 
binaries with M2 = IMq^Mq) and initial distance d = O.Olpc, the mean asymptotic ejection velocity of the primary 
and secondary stars are respectively: v\ ~ 1000km s _1 (1600km s _1 ) and V2 ~ 2800km s _1 (2700km s _1 ). When the 
initial apoapsis of the binary is increased to d = O.lpc we found: v\ ~ 1350km s _1 (2180km s _1 ) and V2 ~ 3500km s _1 
(3200kms _1 ). Table [7l gives the fraction of collisions, mergers, and HVSs (i? e j > lOOOkms -1 ) in the simulations. 
Note that in this table any merger is also counted as a collision. In the table we report the probability of ejection, 
distinguishing between the two components of the binary. It is clear that the initial distance of the binary from the 
central black hole plays a fundamental role in determining which member is ejected. For large initial distances (i.e., 
d = O.lpc) the ejection probability is almost independent on the stellar mass, while for d = O.Olpc, the lighter star is 
preferentially ejected. These results are consistent with the findings of iSari et al.l (2010) that used an approximated 
method to study the dynamical evolution of binaries on parabolic orbits. Our simulations suggest that, in the limit 
that the external orbital energy of the binary goes to zero, the ejection probability becomes an independent function 
of the stellar mass. 
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